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1. INTRODUCTION 

Solar energy is one of the promising types of renewable energy resources in Malaysia with a 
potential total installed solar electricity capacity of more than 6500 MW [1]. The solar electicity is generated 
using either the Grid-Connected Photovoltaic (GCPV) systems or Stand-Alone Photovoltaic (SAPV) 
systems. Mostly, GCPV systems which are widely installed in urban areas where the supply from utility grid 
is highly available. An SAPV system typically consists of PV arrays, solar batteries as well as power 
conditioning units such as charge contollers and inverters. Besides that, the system could also be hybridized 
with other energy resources such as wind, hydro and diesel generator. Moreover, the type, capacity, size and 
configurations of these system components are commonly selected through a systematic process known as 
sizing [1], [2]. For PV-battery-diesel generator normally called as a HSAPV systems. 

HSAPV system usually consists of PV modules, such as power conditioning unit charge controller 
and inverter, a battery bank for storage charges and diesel generator (DG). The size of the HSAPV system 
appears to be more complex and difficult than the size of the photovoltaic system grid-connected since more 
components should be considered in the HSAPV system. In [3], [4], the sizing parameters of the HSAPV 
system is more complicated since the system cannot receive immediate back-up from the utility grid. The 
size of the HSAPV system usually requires system designers to initially choose PV modules, charge 
controller, battery bank, inverter and diesel generator before trying to match the electrical characteristics 
among these components. 
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However, if there is more than one HSAPV model, the sizing process is repeated to achieve the 
most suitable parameters for the design. The conventional sizing process has experienced computational 
burden as reported in [5], [6]. Several studies have been conducted to accelerate the overall sizing process in 
HSAPV design. Most approaches made use of various types of Computational Intelligence (CI) techniques. 
One of the popular optimization techniques that have been widely employed is Artificial Intelligence (AI) 
[7]. AI technique has dealt with difficult problems in some areas; which in turns alleviated the setback 
experienced in the existing conventional techniques. 

There are three important components in AI namely, the Artificial Neural Network (ANN), Fuzzy 
Logic (FL) and Evolutionary Computation (EC) [8]. The three main branches of AI have their own different 
applications.On the other hand, EC is commonly used for population-based metaheuristic optimization 
process. Generally, there are four main branches for Evolutionary Computation (EC) which are Genetic 
Algorithm (GA), Genetic Programming (GP), Evolutionary Strategy (ES) and Evolutionary Programming 
(EP) [9]. In addition, EP and DEA are categorized as the population-based techniques. These meta-heuristic 
methods are also classified as part of CI. 

This paper presents a Hybrid Standalone Photovoltaic System sizing optimization based on load 
profile. There are three optimization techniques have been developed, used to determine the optimum PV 
modules, charge controllers, inverters, battery and diesel generator. Comparative studies have been 
performed FEP, CEP and DEA; while, the ISA used as the benchmark technique. 


2. SYSTEM DESCRIPTION 
This section illustrates the method for sizing HSAPV system. The HSAPV system seen in this 
research is a PV-battery-ganset power system with AC load as shown in Figure 1. 


Charge 
Controller 










AC Load 












Fuel Genset 





Battery Bank 


Figure 1. Block diagram of HSAPV system 


The main thing before forming a HSAPV system, load demand should be determined first. In this 
research, data were collected from a rural area in Kalabakan, Sabah, Malaysia. The estimated daily load 
demand and monthly solar radiations are shown in Table | and Table 2 respectively. 


Table 1. Daily Load Profile 








Anpliance Number of Power Power Usage Bice Usage Ener Surge 
PP units per unit factor Times BY Times By factor 
Dry season Wet season 
(jan..jun) (jul..dec) 
240 ac, SOHz WwW WwW h Wh h Wh 
Water pump 3 120 0.7 2 720 1 360 5 
Ineandeeen 15 40 1 10 6000 13 7800 1 
bulb 
Television 3 80 0.95 8 1920 6 1440 1 
Stand fan 3 60 0.6 4 720 10 1800 5 
Daily Average AC Load Demand (Wh) 9360 11400 
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Table 2. Expected Mothly Solar Irradiation [11] 
month jan feb mar apr may jun jul aug sept oct nov dec 
Irradiation, 
G (Wh/m’) 








4000 4320 4400 4460 4500 4700 5100 5180 5100 5090 5110 5050 





2.1. System Determination 
The steps for sizing the SAPV system are defined as follows [10]: 
Step 1: Determine the required total energy demand, EFveguired_daily 


— Fac 
E+ equired_daily > 7 (1) 


where Eac is the total energy required daily for AC load as delivered by the battery in Wh, 7invis the average 
efficiency of the inverter used for the designed AC loads. 

Step 2: Calculate the System Voltage (SV) for the HSAPV system. The SV can be selected as 12V, 24V or 
48V depending on size of load demand required as shown in Table 3. 


Table 3. Range of SV [10] 











tes Recommended 
Condition SV (V) 
Case 1 Eequired daily < 1kWh 12 
Case 2 1kWh < Ejequired_daily S 4 kWh 24 
Case 3 Erequired daily = 4 kWh 48 
Step 3: Determine the total capacity required daily, Crequired_daily 
C m . a Exequired_daily (2) 
required_daily — SV 
Step 4: Determine the battery bank capacity required, Chank_required 
C =C Tautonomy 3 
bank_required — ‘required_daily x DODmax ( ) 


where DOD ax is the maximum battery depth of discharge. Tautonomy represents the expected number of days 
that the load demand will be happened by battery bank without sunlight availability. 

Step 5: Determine the total load current from the battery bank, [hotai_toad_curren and the battery bank discharge 
rate, Thank disch 


¥ ACPower 


(4) 


1 
Total load_current ety Dx DCPower + 
sv Tiny 


Chank_required (5) 


Thank_disch = I 
total_load_current 


where Thank isch iS the discharge rate of battery bank in h. 
Step 6: Determine the battery bank configuration, i.e. the number of battery in series string, Nseries bank, 
number of battery strings in parallel, Nparaitel_bank and bank capacity of the selected battery, Coank_setected, 


SV 


Nseries _bank “7 (6) 
nom_batt 


where Vnom_batt is the nominal battery voltage and Cper battery is the AH capacity which can selected from the 
datasheet 


Crevised_bank_required 
N aa—arraearrr (7) 


llel_bank — 
DOPE. aR, Cyer_battery 


Cpank_selected = parallel_bank x Cyer_battery (8) 
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Step 7: Based on MPPT charge controller; determine the maximum and minimum number of series PV 
module based on open circuit voltage. Nseries_max_based_on_Voc and Neries_min_based_on_Voc. 


N _ 0.95 X Vmax cc 
series_max _based_on_Voc — 


(9) 


Vmax _oc 


24 *V min _win.ce (10) 


Nseries_min_based_on.Voc = Vai 
min _mp x fpov_ce 


where Vinax_cc is the maximum input voltage rating of the charge controller and Vinax_oc is the maximum open 
circuit voltage rating of PV module. Vinin_window_cc is the minimum allowable window input voltage to the 
MPPT of the charge controller, Vinin_mp is the minimum voltage at maximum power of PV module and fpy cc is 
cabling efficiency from PV to charge controller. 

Step 8: Determine the maximum number of series PV module based on maximum power voltage, 


N, max_series_based_on_Vmp 


0.95 x Vmax win_cc (11) 


Nseries _max_based_on_Vmp = V, 
max _mp 


where Vimax_win_cc 18 the maximum window voltage of the charge controller and Vinaxmp is the maximum 


voltage at maximum power at PV module. 
Step 9: Determine the corrected output power of module, Prod corrected and sub-system efficiency, Pmod_corrected 


Pmod_corrected = Pmp_ste x hemp x fm x fairt (12) 


where Pymp sic is the rated power at STC, fair: is the derating factor for dirt, finn is the derating factor of the 
manufacturing tolerance and fiemp is the derating factor for temperature. 


Iyvy_ss ~ (13) 


I nv_cable x 1 controller x Myatt 


where 7pv cable is the efficiency of cable from PV array to battery, 7controler is the efficiency of the MPPT 
charge controller and 7pat is the efficiency of the battery. 
Step 10: Determine the total number of PV modules for the PV array, N7_mpp: 


Eyequired_daily* fo 
Nop omppt = (14) 
= Pmod_corrected X PSH x Tyv_ss 


where 7pv ss is the sub-system efficiency of PV array and Pip sic is the rated maximum power at Standard 
Test Conditions (STC). 
Step 11: Determine the number of parallel PV modules, Np_npp:, 


NT_mppt (15) 


N = 
P_mppt Ns mppt 


where Ns mppr is the total number of PV modules in series between Neeries_max ANd Neeries_min 
Step 12: Determine the peak capacity of array using the MPPT charge controller, Parray_mppt_ce 


Parray_mppt.cc = Pmp mod_ste x Nz mppt x Np mppt (16) 
Step 13: Determine the apparent power of the inverter required for 30 minutes at maximum demand. 


Sinv_30min = Ymax_AC_demand X Sfiny (17) 


Step 14: Determine the apparent power of the inverter required for surge demand. 





Sinv_surge = Smax_AC_surge x Sfinv (18) 
Prea 
Smax_AC_surge = ma X Sy (19) 


where Siny_3omin iS the apparent power of the inverter at maximum demand, Sin surge is the surge rating of the 
inverter, Simax_AC_surge is AC the surge load demand, Smax_Ac_demandis the maximum AC load demand, Sfiny is the 
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safety factor of inverter, P,eaj is the real power consumption when running, S; is the surge factor and PF is the 
power factor. 
Step 15: Determine the minimum size of diesel generator; 


— (Smax_AC_demand~Sinv_30min) 
Sgen_max_demand a facrare X figo (20) 
where Syen_max_demand_parallel is the minimum diesel generator rating based on power demand, S max_Ac demand 18 
the maximum AC load, Sinv_3omin is the 30 minutes’ apparent power rating of the inverter, f,. is diesel 
generator factor and fierate is the total derating factor. Sgen_max_surge can be calculated using Equation (21). 


oe (Smax_AC_surge—Sinv_surge) 
Sgen max_surge —~ X figo (21) 
5 - fderate 


where S'gen_max_surge 18 the minimum diesel generator rating based on surge condition, Smax_Ac_surge is the AC 
surge load demand and Siny_suree is the surge rating of the inverter. 


(Sbatt_charger—Smax_AC_demand) 
Sgen_batt a X fo (22) 


fderate 
where Syen_batt 18 the diesel generator rating based on battery charging and maximum load profile. 
Step 16: Calculate the Performance Ratio, PR 


PR = ——25eee —_ (23) 


Parray_stc x PSHanual 


where E’sys exp is the expected kWh output of the system and Parray sicis the nominal PV array power at STC. 
PSH represents the estimated annual PSH derived from Table 2. 


3. DOLPHIN ECHOLOCATION ALGORITHM 

Natural resources have inspired many of the man-made technologies. Researchers believe that 
dolphins are the second to humans in smartness and intelligence. Dolphin can use the echolocation method to 
help them navigate better in underwater such as to find foods or communicating. The sounds which produced 
from a dolphin are in the form of clicks which has higher frequency than the sound for communication. 
These clicks create sound waves that able to travel quickly and when the sound waves bounce off objects 
around them, they return to the dolphins as echoes and the dolphins can determine the distance [12]. 

Basically, the DEA-based optimization process, it starts with the initialization of an objective 
function. After that, the Predefined Probability of a loop, PP (Loop;) and the Accumulative Fitness of the iy, 
alternative to be chosen for the j,, variable, AFjj are calculated using 


PP(Loop;) = PP(Loop,) + (1 — PP(Loop,)) x (or) (24) 
te 
Re 


AF;; = x (R, — |k|) x Fitness @) + AFij-1 (25) 
where PP(Loop;) is the convergence factor of the first loop. Loop; is the number of current loop and 
LoopNumber is the number of loops. AFjj.; is the previous AF of the i, alternative to be chosen for the j,;, 
variable. Besides, R, is the valuable radius where the AF of the i,, alternative to be chosen for the j,, variable 
is affected from its fitness. This radius is recommended to be not more than quarter of the search space [13]. 
k is the positive and negative value for the effective radius where the AF of the i, alternative to be chosen for 
the jin variable are affected from its fitness and Fitness (i) is the fitness of the location i. 

Next, a small value of ¢ is added to AF are consistently in the search space for the selected option. 
Up to now, previous studies have indicated that, the value of € is important and superior to be less than the 
minimum value achieved for the fitness [13]. 


AF =AF +e (26) 


After that, set the AF=0 for the best candidates from the population while the rest of the candidates 
are maintained with the original value for their AF’. The step can be defined as follows: 
For j=1: Number of variables ; i=1: Number of alternatives; If i=The best location (j) 


AF, = 0 (27) 
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Next step is to determine the probability, P;; of choosing the i alternative for all the number of 
variables, j is calculated using: 


Pij = ray (23) 


After that, PP(Loop;) is allocated to all alternatives chosen for all variables. Then, the next location 
of candidate is determined according to the probabilities assigned to each alternative. In this study, the PR 
value was determined based on the fitness value of each candidate. The results were ranked accordingly to 
their fitness values with the candidate having maximum PR. Finally, all steps are repeated until the maximum 
loop numbers is reached. 














Selection of components HSAPV 
system 
Sizing of HSAPV system 


Installation of the candidates 
Calculate P(loopi) for every loop 
Accumulative fitness (AF) 
calculation 
Add all AF with ¢ 
Find the current best 


$ the populatio 
tthe current best 












Maintain the AF 


Calculate the probability of 
choosing AF 





Allocate the probability of 
choosing AF 


Calculate the next step locations 
Fitness calculation 
Rank and select the highest fitness 
value 







Figure 2. Flowchart of sizing and optimization technique using DEA 
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4. RESULTS AND ANALYSIS 

This section describes the results of HSAPV system using DEA and EP models with ISA as shown 
in Table 4. These results finding the maximized PR and the optimal configuration among set of system 
components which can meet the load demand requirements. From the table, it is observed that all the four 
technique ISA, DEA, CEP and FEP conclude no of inverter code is 7. Beside that, the number of charge 
controller code identified by ISA, DEA and FEP is identical, i.e. 7; while, CEP determines 8 as the charge 
controller code value. ISA is the technique for the benchmarking of the whole optimization algorithm 
techniques. Apart from that, the performance DEA and FEP was found to produce the optimal value of PR 
with ISA, i.e. 0.7682. However, DEA has speed up the computation time to accomplish the minimal value. 
This result shows 91.32% faster than ISA technique. 


Table 4. Performance of Comparison of DEA and EP Models with ISA 








EP Techniques 
Results ISA CEP FEP DEA 
Battery code 12 12 12 12 
PV code 12 6 12 12 
Inverter code 7 7 7 7 
Charge Controller code 7 8 7 7 
Diesel Generator Code 11 11 11 11 
N_batt_s 8 8 8 8 
N_batt_p 2 2 2 2 
N_pv_s 2 3 2 2 
N_pv_p 8 6 8 8 
Nt_inv 1 1 1 1 
Nt_cc 2, 2 2 2 
Nt_diesel_gen 1 1 1 1 
Optimal, PR 0.7682 0.7428 0.7682 0.7682 
Comp. time, (s) 2003.04 174.58 174.02 173.91 





The next analysis is the maximization of PR based on four techniques by using different population 
size is tabulated in Figure 3. Based on figure, for the small size population (n=5), the range of maximum PR 
for three optimization is 0.71-0.750 compare with ISA (i.e 0.7628). Furthermore, when the population size is 
increased, the values of fitness values are increase accordingly. On the other hand, the FEP and DEA reached 
the maximum value of PR (i.e 0.7628) with population size of 30 as those set by ISA, while the CEP is 
reached at value 0.7448 until the end of population size. In addition, the DEA is maintained with the 
maximum value of PR at population size of 25 onwards (i.e 0.7628), while the FEP is reached the maximum 
PR at population size of 30. These mean that DEA is better than CEP and FEP in optimizing the HSAPV 
system for maximization of fitness value. 
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Figure 3. Performance of fitness value with different techniques at different number of population 


5. CONCLUSION 

This paper has presented Hybrid Standalone Photovoltaic Systems Sizing Optimization based on 
load profile. There are three different optimization techniques were investigated before selecting the best 
optimal value based on sizing algorithm. From the results, the DEA has reached the maximum of fitness 
values, together with the minimal computation time. The DEA has been compared with the ISA, whereas 
ISA is the iterative technique where can produced the optimal value of PR. The development of DEA is also 
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feasible for solving the other optimization problems. For future works, DEA could be potentially conducted 
to apply in multisource of SAPV systems. 
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